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We have studied the diamond— > /3-tin phase transition in Si using diffusion quantum Monte Carlo 
(DMC) methods. Slater-Jastrow-backflow trial wave functions give lower DMC energies than Slater- 
Jastrow ones, and backflow slightly favors the /3-tin phase with respect to the diamond phase. We 
have investigated the changes in the equation of state that result from the use of different pseu- 
dopotentials, the inclusion of either zero-point motion or finite-temperature vibrations, and the 
application of corrections for finite-size effects. Our tests indicate that the choice of pseudopotential 
can significantly affect the equation of state. Using a Dirac-Fock pseudopotential leads to an over- 
estimation of the transition pressure but an empirical pseudopotential designed for use in correlated 
calculations gives a transition pressure in quite good agreement with experiment. 

PACS numbers: 02.70.Ss,71.15.Nc,64.60.Ej 

I. INTRODUCTION 

Calculating accurate energy differences between the atomic structures of materials is one of the central problems 
in quantum mechanical studies of condensed matter. Here we study the first-order pressure-induced structural phase 
transition between the diamond and /3-tin phases of Si. This is a demanding problem because of the substantial volume 
change and the increase in the coordination number from four to six at the transition, and the fact that the structural 
phase transition is accompanied by a metal-insulator transition. Many theoretical studies of the diamond— 5- /3-tin 
phase transition in Si have been reported, mainly using first-principles density-functional theory (DFT) methods*^— 
The DFT results show significant dependence on the density functional used^ indicating the necessity for an accurate 
description of exchange-correlation effects. A similar dependence on the functional has been noted for self-interstitial 
defects in Si3~— These results motivate studies using other descriptions of many-body effects such as that provided 
by quantum Monte Carlo (QMC) methods. 

Pressure-induced phase transitions have been studied using the diffusion quantum Monte Carlo (DMC) 12 ' 13 and 
auxiliary-field quantum Monte Carlo (AFQMC) methods.— The pressure-induced phase transition of Si from the 
diamond to the /3-tin structure has been studied previously using both DM C 15 ' 16 and AFQMC ^ while DMC studies 
of structural phase transitions have also been reported for N2r^ MgOp^ FeO^ Si02rl and BNi^ 2 - 

The purpose of the present study is to test recent developments in DMC methodology We have used Slater- 
Jastrow-backflow (SJB) trial wave functions which go beyond the single-particle nodal surfaces of Slater- Jastrow (SJ) 
wave functions and significantly lower the energies of the diamond and /3-tin phases. We have tested two recent and 
very different schemes for calculating corrections to the DMC energies obtained with finite supercells that allow us 
to extrapolate to the infinite-cell limit. To compare the DMC transition pressure and volumes with experiment it is 
necessary to include the vibrational free energy, which wc have re-calculated within DFT with greater accuracy than 
was achieved in previous studies. Finally we have calculated the energy-volume relations within DFT using different 
pseudopotcntials in order to study the likely influence of the pseudopotential on the DMC transition pressure and 
volumes and the equations of state of the diamond and /3-tin structures. 

The rest of the paper is organized as follows. In Sec.[II]we define the Hamiltonians used and in Sec. lIIII the trial wave 
functions are described. Finite-size corrections are discussed in Sec. IIVI and the phonon calculations are described 
in Sec. [Vj The main DMC results are reported in Sec. I VII Section IVIII describes our method for testing different 
pseudopotcntials and reports and discusses results for the cquation-of-state parameters and transition pressures and 
volumes. We draw our conclusions in Sec. IVIIII The QMC calculations were performed with the CASINO code^ 3 - and 
for the DFT calculations we used the CASTEP plane-wave pseudopotential code^ 
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II. THE HAMILTONIAN 



We used the Dirac-Fock pseudopotential of Trail and Needs^ 5 - to describe the interaction of the valence electrons 
with the Si 4+ ions. We also added a core-polarization potential (CPP) to each ion i 26 ' 27 In the CPP approximation, 
the polarization of each core is determined by the electric field at the nucleus from the instantaneous positions of 
the electrons and the other atomic cores. CPPs therefore give an approximate description of both dynamical core- 
valence correlation and static polarization effects. We used the CPP parameters from Ref. [2fj| and the CPP energy 
was evaluated using the scheme described in Ref. l28l The non-local pseudopotential energy was calculated using the 
variational scheme of Ref. I29I . 

We used 54-atom simulation cells subject to periodic boundary conditions. To make the simulation cell for the 
diamond structure we took the standard set of primitive translation vectors of the two-atom face-centered-cubic cell 
and multiplied them by three to give a cell accommodating 54 atoms. The /3-tin structure has a body-centered 
tetragonal lattice and a two-atom primitive unit cell. We used a c/a ratio of 0.556 for the QMC calculations, which 
is close to the experimental value of 0.5498 reported at the phase transition^ The energy difference between the 
structures when using these two c/a ratios is negligible on the scale of interest. The primitive lattice vectors of the 
two-atom /3-tin structure can be chosen to be u = (a, a, c), v = (a, —a, c), and w = (a, a, — c). We wrote a code to 
search over all possible 54-atom cells constructed from linear combinations of u, v, and w and chose the cell with the 
maximum value of the shortest distance between supercell lattice points, which has translation vectors U, V, and W, 
where U = — 9u + 3v + 6w, V = 3v, and W = — 3u + 2v + w. Using such a simulation cell gives the largest separation 
between each electron and its periodic image, which helps to reduce the finite-size effects arising from long-ranged 
exchange-correlation effects and allows the isotropic terms in the Jastrow factor and backflow transformation (see Sec. 
IIIIP to give a good description of the electron correlation. The Ewald potential was used to model the interactions 
between point charges. 



III. TRIAL WAVE FUNCTIONS 



The DMC method relies on the availability of an accurate trial wave function which is used to guide the sampling 
and to enforce the fixed-node constraint ! 12 ' 13 We used trial wave functions of SJ and SJB types. The Slater- Jastrow 
form is 



*sj(R)= cxp[J(R)] det Vn(rJ) 
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where R denotes the positions of all the electrons, rj is the position of the ith spin-up electron, rj is the position of 
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are determinants of 



the jth spin-down electron, exp[J(R)] is the Jastrow factor, and det ?/> n (rJ) 
single-particle orbitals. 

The orbitals were calculated using the CASTEP code and a large basis-set energy cutoff of 1088 eV. The orbitals 
were then transformed into a B-spline or "blip" polynomial basis for greater efficiency^ We used Jastrow factors 
consisting of polynomial electron-nucleus (en) and electron-electron (ee) terms and a plane-wave ee term^ 2 - with a 
total of 26 optimizable parameters. We used wave functions both with and without backflow transformations^ In 
the SJB wave function the particle coordinates in the determinants in Eq. (TTJ) are replaced by collective coordinates 



x»(R) 
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where £j(R) is the backflow displacement of electron i, which depends on the positions of all of the electrons and 
ions. The backflow displacement consisted of polynomial ee and en terms^ 3 . and contained a total of 24 optimizable 
parameters. The wave-function parameters were optimized by a variational Monte Carlo (VMC) procedure in which 
we first minimized the variance of the energ y 34 ' 35 and then minimized the variational energy itself^ The wave function 
parameters were optimized at a single k s point in the Brillouin zone of the simulation cell (this being L and M for 
the diamond and /3-tin structures, respectively) and the same Jastrow and backflow parameters were subsequently 
used for each k s point studied, where k s is the many-body Bloch vector^ 

All of the results reported here were obtained using a DMC time step At of 0.01 a.u. Some test calculations using a 
SJ trial wave function were also performed with At =0.02 a.u., which give very similar results to those with At =0.01 
a.u. We used a target population of 640 configurations in each of our DMC calculations. 



3 



0.14 




^0.08 

Q 

hf 0.07 



100 120 140 160 
Volume (a.u. / atom) 



200 



FIG. 1: (Color online) Difference between the SJ and SJB DMC energies for the diamond and /3-tin phases. The diamond- and 
/?-tin-structure calculations were performed at the L and M points in the Brillouin zone, respectively. 



IV. FINITE-SIZE CORRECTIONS 



The 54-atom simulations give a reasonable description of the infinite system, but it is still necessary to make finite- 
size corrections to obtain accurate results. The finite-size effects can be divided into a part which is very similar 
to the Brillouin zone integrations performed in single-particle calculations, and exchange-correlation effects arising 
from electron separations beyond the simulation cell. Brillouin zone integration is more important in metals, where 
it is required to describe Fermi-surface effects, than in insulators, where a sparse sampling of the Brillouin zone 
is sufficient. Our 54-atom cells are equivalent to sampling the Brillouin zone of the two-atom primitive cell in a 
single-particle calculation with 27 k points. 

It is only possible to perform calculations with one simulation-cell Bloch vector k s at a time within explicitly 
correlated methods. For the insulating diamond phase we used only the L-point, which allows the use of a real trial 
wave function while giving a good description of the energy and charge density in single-particle calculations! 37 ' 38 
For the metallic /3-tin phase we performed DMC calculations at twelve randomly chosen k s points and averaged the 
results^ 

We tested two different schemes for calculating additional finite-size corrections. In the first scheme we added 
a further k-point correction consisting of the difference between the PBE-DFT energies with a very dense k-point 
sampling and with the finite k-point sampling corresponding to that used in the QMC calculations. We then added 
a correction arising from the long-ranged exchange-correlation interaction using the method described in Refs. 
and l4ll . which involves evaluating the static structure factor and two-body Jastrow factor on the grid of reciprocal 
lattice vectors within QMC and interpolating to obtain an approximate static structure factor and two-body Jastrow 
factor for the infinite system, from which a finite-size correction can be evaluated. The second scheme is that 
of Kwee et a/42 which involves performing a DFT calculation using a modified local density approximation (LDA) 
functional obtained from DMC calculations for homogeneous electron gases in finite simulation cells subject to periodic 
boundary conditions. This method includes both the k-point and long-ranged exchange-correlation corrections, but 
it also involves making the somewhat stronger assumptions that the LDA is sufficiently accurate and the correction 
is insensitive to the shape of the simulation cell. 

We found excellent agreement between the finite-size corrections obtained from the structure factor metho d 40 ' 41 
and the DFT method of Kwee et al^ The finit e-size corrections themselves are quite large, amounting to an increase 
in energy of about 0.11 eV per atom for the diamond structure and about 0.099 eV per atom for /3-tin. However, the 
difference between the finite-size corrections obtained from the two schemes is, in the worst case, about 0.007 eV per 
atom, while the average of the absolute deviation over all of the data points for the two structures is 0.004 eV per 
atom. 



V. PHONON CALCULATIONS 



We evaluated the phonon contribution to the Hclmholtz free energy of each phase as a function of volume using the 
PBE functional. We used a 64-atom supercell, an energy cutoff of 653 eV, and a 6 x 6 x 6 k-point mesh to calculate 
the phonons of the /3-tin structure and a 54-atom supercell, a 435 eV cutoff, and a 3 x 3 x 3 k-point grid for the 
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FIG. 2: (Color online) Variation of the free energy with volume for the diamond and /3-tin phases of Si calculated within DMC 
and with corrections for the vibrational energy from PBE-DFT. Static lattice, T — and T = 300 K results are shown, and 
the zero of free energy is taken to be the DMC results including vibrational effects at 300 K. The common tangents labeled 
"Tang." join the coexistence points. Filled symbols denote that the structures were found to be dynamically stable within the 
DFT phonon calculations. The statistical error bars on the DMC energies are in the range 0.00005-0.00007 a.u. per atom in 
all cases. 

diamond structure. The phonons frequencies and modes were calculated using a finite-displacement method and the 
results were extrapolated to zero displacement. Gaal-Nagy et al£2- have shown that it is difficult to obtain stable 
phonon modes of the /3-tin phase of Si in DFT calculations. Our study supports this conclusion, and we have found 
stable phonon modes only by carefully relaxing the c/a ratio at every volume and, even then, the phonon modes are 
only stable over the volume range 91.4-111.8 a.u. per atom. Over the volume range where the phonons are stable, the 
variation of the relaxed c/a ratio with volume is well fitted by c/a = 0.521075 + 0.377765/(131.41 - y) 77666 , where 
V is the volume per atom in a.u. 

The phonon free energy was calculated from the phonon frequencies. The previous DMC and AFQMC calculations 
estimated the changes in the transition pressure from nuclear zero-point motion and finite-temperature vibrations 
using the DFT vibrational data of Gaal-Nagy et al~ These data suggest that including zero-point motion reduces 
the phase transition pressure by about 1 GPa, and that finite-temperature effects at 300 K reduce the transition 
pressure by a further 0.3 GPa. We believe that Gaal-Nagy et al£^ have somewhat overestimated the size of the 
vibrational effects because they used a non-optimal c/a ratio. Our phonon calculations give a reduction in the 
transition pressure from zero point motion of about 0.62 GPa and a further reduction of 0.34 GPa at 300 K. These 
values do not differ significantly when we use LDA, PBE, or DMC static-lattice energy-volume curves. 



The reduction in the DMC energy resulting from the inclusion of backflow is shown in Fig. Q] The energy reduction 
is greater in the /3-tin phase than in the diamond phase at the same volume, which we conjecture is due to the metallic 
nature of the /3-tin phase. The energy reduction increases with the density of the system, as has been found in studies 
of the electron gas4£ The use of backflow reduces the energies of the diamond structure by 0.07-0.1 eV per atom and 
the /3-tin structure by 0.11-0.13 eV per atom, over the respective volume ranges, but the energy difference between 
the structures at the transition is reduced by only about 0.015 eV. 

The total Helmholtz free energy at temperature T at each volume was evaluated as the sum of the DMC energy, 
the finite-size correction from the scheme of Kwee et al.^ and the phonon free energy. The parameters of the Vinct 
equation of stated were then calculated using a least-squares fitting procedure. The Helmholtz free energies for the 
static lattice, zero temperature, and T = 300 K cases are shown in Fig. [2j The fitted equation-of-state parameters 
for the diamond and /3-tin structures are given in Tables HI and HT1 respectively, and coexistence pressures, etc., for the 
phase transition are reported in Table IIIII 



VI. DMC RESULTS 
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Source Vo 


(a.u. / atom) 


B (GPa) 


B' 


Experiment 


135. l a 


97.9(1) 6 


4.24 b 


DMC 


134.3(1) 


97.1(3) 


4.18(5) 


DMC + BFD-pp 


132.3(1) 


99.7(3) 


4.18(5) 


DMC + EMP-pp 


133.3(1) 


96.2(3) 


4.19(5) 


DMC + PBE-pp 


153.7(1) 


85.3(6) 


3.67(5) 


DMC + LDA-pp 


131.4(1) 


101.5(3) 


4.18(5) 


PBE-pp 


138.0 


86.0 


4.52 


LDA-pp 


132.8 


94.4 


4.46 



°Ref. [H 
b ReL [H. 



TABLE I: The equilibrium volume per atom Vo, bulk modulus Bo, and pressure derivative of the bulk modulus B' of diamond- 
structure Si from experiment and from DMC and DFT calculations. All results are for 300 K. 



Source 


V (a.u. / atom) B (GPa) 


B'o 


DMC 


103.1(1) 


102.7(26) 


4.7(2) 


DMC + BFD-pp 


101.4(1) 


106.2(26) 


4.7(2) 


DMC + EMP-pp 


102.4(1) 


104.2(26) 


4.7(2) 


DMC + PBE-pp 


101.5(1) 


105.3(21) 


4.7(2) 


DMC + LDA-pp 


100.6(1) 


107.4(21) 


4.8(2) 


PBE-pp 


104.8 


94.1 


4.74 


LDA-pp 


103.9 


96.0 


4.76 



TABLE II: The equilibrium volume per atom Vo, bulk modulus Bo, and pressure derivative of the bulk modulus B' of /3-tin- 
structurc Si from DMC and DFT calculations. All results are for 300 K. 



Source pt (GPa) Vdia (a.u. / atom) V^g-tm (a.u. / atom) AVt (a.u. / atom) AFo (eV / atom) 



Experiment (Ref. 30) 




11.7 


122.48 


94.17 


-28.31 




LDA (This work) 




6.17 


125.42 


95.63 


-29.79 


0.2681 


PBE-GGA (This work) 




8.75 


126.82 


96.95 


-29.87 


0.2538 


DMC (This work) 




15.1(1) 


119.0(1) 


92.1(1) 


-26.9(1) 


0.405(3) 


DMC + BFD-pp (This 


work) 


15.7(1) 


117.1(1) 


90.5(1) 


-26.6(1) 


0.418(3) 


DMC + EMP-pp (This 


work) 


12.2(1) 


120.3(1) 


93.2(1) 


-27.1(1) 


0.329(3) 


DMC + PBE-pp (This 


work) 


14.9(1) 


134.1(1) 


91.0(1) 


-43.1(1) 


0.656(3) 


DMC + LDA-pp (This 


work) 


16.9(1) 


115.6(1) 


89.3(1) 


-26.3(1) 


0.446(3) 


DMC (Ref. ^5) 




15.5(5) 










DMC (Ref. 16) 




14(1) 


122.4(3) 


93.8(7) 


-28.69(8) 


0.42(2) f 


AFQMC (Ref. _17) 




12.6(3) 


122.48* 


94.17* 


-28.31* 





TABLE III: The experimental diamond— >• /3-tin phase transition pressure and volumes and calculated coexistence pressures and 
volumes in Si, and the difference between the zero-pressure free energies. Volumes are given in a.u. per atom. The asterisks 
denote that the experimental value was assumed. The data are for T = 300 K except for the value marked with a dagger, 
which is for zero temperature. The addition of a 300 K temperature correction to the value of AFo from Ref. [l6| would result 
in a decrease of about 0.024 eV per atom. 



VII. TESTING DIFFERENT PSEUDOPOTENTIALS 



Pscudopotcntials are very well established in DFT calculations, where they can give results in excellent agreement 
with all-electron calculations. The accuracy of pseudopotentials in reproducing all-electron results in explicitly corre- 
lated calculations such as QMC is rather more doubtful. One difficulty is that there is no straightforward equivalent 
of the method normally used within independent electron theories of constructing a pseudopotential by inverting 
the atomic Schrodingcr equation for the valence electrons. Indeed, the idea of dividing the electrons into core and 
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valence shells is not well defined in the many-body system. The difficulties in generating pscudopotcntials which are 
consistent with correlated methods for calculating the valence electronic structure leads to additional errors. 

To examine the DMC results using different pscudopotcntials we calculated corrections to the energy-volume curves 
using different pseudopotentials. We modeled the interactions between the valence electrons by the PBE-GGA func- 
tional and approximated the DMC energy for pseudopotential X-pp by 

pX-pp _ pDF-pp , / pX-pp _ pDF-pp\ / o \ 

DMC — DMC ' ^ PBE PBE J • \°> 

where the subscript refers to the method used to calculate the energy, and DF-pp denotes the Dirac-Fock 
pseudopotential^ We tested four pseudopotentials: ultrasoft LDA and PBE pseudopotentials, the empirical (EMP) 
pseudopotential of Ref. [4g|, and the Hartree-Fock pseudopotential of Burkatzki, Filippi, and Dolg (BFD) from Ref. 

There is some evidence that Hartree-Fock (and presumably Dirac-Fock) pseudopotentials give better results within 
DMC than DFT pseudopotcntials^i The EMP pseudopotential uses the measured values of the ionization energy of 
Si +3 ions which have a single electron in the 3s, 3p, or 3d orbital outside of the neon core, and calculated data for 
the amount of valence charge density outside of the pseudopotential core radius in the ionic states. Corrections for 
the core relaxation effects arising from generating the pseudopotentials in ionized configurations arc also included. 
Note that Esler et al£l have developed a method for calculating a correction to the QMC energy that accounts for 
the error due to the pseudopotential approximation. The method uses all-electron results for small systems, although 
these would be very costly to evaluate for Si. 

A. Comparison of equation-of-state parameters 

Data for the equilibrium volume Vo> the bulk modulus Bq, and the pressure derivative of the bulk modulus B' of 
diamond-structure Si at 300 K are given in Table HI As shown in numerous DFT calculations, the equilibrium volume 
Vq of diamond-structure Si with the PBE functional is larger than experiment while the LDA value is smaller than 
experiment. The PBE value of Bq is significantly smaller than experiment while the LDA value is close to experiment, 
and the PBE and LDA values of B' are rather larger than experiment. 47 ' 48 The experimental value of Bq is believed to 
be very reliable and therefore we can be confident that PBE gives too small a value. Reference |48| does not report an 
error bar for the measurement of B' and so it may not be safe to draw conclusions about the accuracy of experimental 
or calculated values for it . 

The DMC calculations for the diamond structure with the Dirac-Fock pseudopotential provide very accurate results 
for Vo and Bq. The DMC calculations with corrections for using other pseudopotentials [see Eq. ([3])] give results 
of variable quality. The equation-of-state parameters obtained using the BFD pseudopotential are very similar to 
those obtained using the Dirac-Fock pseudopotential. This is to be expected because these pseudopotentials were 
designed to reproduce Hartree-Fock or Dirac-Fock results, and the rclativistic effects included in the Dirac-Fock 
pseudopotential are unimportant for Si. However, the methods used to construct the pscudopotcntials are rather 
different, and it is therefore relevant to compare the results obtained with them. An earlier comparison of the Trail- 
Needs Dirac-Fock pseudopotentials and the BFD pseudopotentials for the LiH molecule also found that very similar 
results are obtained with the two classes of pseudopotential^ The equation-of-state parameters obtained with the 
EMP pseudopotential are close to those from the Dirac-Fock pseudopotential, but the results obtained with the PBE 
and LDA pseudopotentials are significantly different and are in poorer agreement with experiment. We note that the 
DMC calculations with the Dirac-Fock, BFD, and EMP pseudopotentials, which give the best values of Vq and Bq, 
also provide the values of B' Q in the best agreement with experiment (note, however our remark about the lack of an 
error bar on the experimental value of B' Q ). The PBE and LDA pscudopotcntials give results which differ substantially 
from experiment. 

The equation-of-state parameters for ^-tin-structure Si at 300 K are given in Table HT1 Experimental data for Vq, 
Bq, and B' are not available for the /3-tin structure of Si. We note that the values of Vq and B' Q for the /3-tin structure 
are similar for the DMC- based and DFT-based calculations, but the DMC-based calculations give larger values of Bq. 

B. Phase transition pressures and volumes 

The best experimental value for the room temperature phase transition pressure is probably that of McMahon and 
Nelmes of pt = 11.7 GPap£ which is close to the value of p t = 11.3(2) GPa found earlier by Hu et al£& Reference l30l 
also reported the structural parameters of the diamond and /3-tin structures at the transition. Several factors should 
be borne in mind when comparing the experimental and theoretical results. The experimental transition pressure 
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and volumes were measured on compression while theoretical studies normally report the coexistence pressure, which 
is expected to be lower because strongly first-order transitions are likely to involve passing over kinetic barriers. 
Unfortunately the diamond— > /3-tin phase transition in Si is irreversible and different structural phases are found 
on decompression^ so that the coexistence pressure cannot be bracketed by the values obtained on compression and 
decompression. Prior to compression, the powdered sample is isotropic, but the compression has a uniaxial component 
and, after the transition, the /3-tin phase shows a "preferred orientation" , that is, the microcrystallitcs have a tendency 
to align, indicating that the stress on the sample during the transition is anisotropic. Various calculations^— have 
shown that an anisotropic stress lowers the transition pressure (the pressure is equal to minus one third of the trace 
of the stress tensor) . As the effects of the anisotropic stress and the kinetic barrier act in opposite directions on the 
measured transition pressure, and it is difficult to estimate their size, it seems reasonable to neglect them. Finally, the 
experimental data were measured at room temperature, while the electronic structure calculations were performed at 
zero temperature and the vibrational effects were added afterwards. 

The data in Table Hill show that the DMC calculations with finite-size and temperature corrections give a coexistence 
pressure of 15.1(1) GPa, which is 3.4 GPa larger than the measured transition pressure and the DMC calculations 
underestimate the transition volumes. Correcting the DMC energies to give an approximation to the results that 
would be obtained using the EMP potential using Eq. ([3]) gives results for the coexistence pressure and volumes in 
quite good agreement with the experimental data. Correcting the DMC energies to approximate the results that 
would be obtained for the PBE and LDA pseudopotentials leads to poorer transition volumes and pressures. 

VIII. CONCLUSIONS 

Backflow improves the description of exchange-correlation effects and significantly lowers the VMC and DMC 
energies of the diamond and /3-tin phases of Si. These energy reductions, however, mostly cancel in the DMC energy 
differences between the two phases and, overall, backflow favors the /3-tin phase with respect to the diamond phase 
by roughly 0.015 eV per atom at the transition. In a previous DMC study it was suggested that the overestimation 
of the experimental transition pressure was most likely to be due to fixed- node errors^ The current much-refined 
study includes backflow transformations which reduce the fixed-node error, but the DMC results with the Dirac-Fock 
pseudopotential (corrected for finite-size and vibrational effects) are within error bars of the results of Ref. [l5|. The 
significant reductions in the energies of the two phases due to the use of backflow make it less likely that the residual 
errors arise from the fixed-node approximation. 

We tested two different schemes designed to correct for the finite sizes of the simulation cells. The finite-size 
corrections arc of order 0.1 cV per atom, but the maximum difference between the corrections provided by the two 
schemes is only about 0.007 eV per atom. This suggests that both schemes are working extremely well in this system. 
We calculated the phonon frequencies of the diamond and /3-tin phases as a function of volume, finding that it is 
very important to relax the c/a ratio of the /3-tin phase at each volume to obtain stable phonon modes. Even then 
the /3-tin phase was found to be dynamically stable only over a small volume range. Our phonon calculations give 
a reduction in the transition pressure from including finitc-tcmpcraturc effects at 300 K of nearly 1 GPa, compared 
with the reduction of 1.3 GPa obtained by Gaal-Nagy et al£^ 

Finally, we estimated the DMC results that would be obtained with different pseudopotentials and found a sub- 
stantial dependence on the pseudopotential. The Dirac-Fock pseudopotential used in our QMC calculations gave 
a transition pressure of 15.1(1) GPa, while the BFD Hartree-Fock pseudopotential gave very similar results. The 
EMP pseudopotential gave 12.2(1) GPa, a PBE pseudopotential gave 14.9(1) GPa, and an LDA pseudopotential gave 
16.9(1) GPa. The best agreement with experiment was obtained with the EMP pseudopotential. It is worth noting 
that the EMP pseudopotential is the only one designed for use in explicitly correlated calculations, and that it gave the 
best results among the Si pseudopotentials tested in correlated calculations in Ref. |49|. The values of the equilibrium 
volume, bulk modulus, and the pressure derivative of the bulk modulus obtained with the EMP pseudopotential were 
also in good agreement with experiment. 

The results of our QMC calculations for Si and, presumably, calculations with other atoms, depend significantly on 
the pseudopotentials used. We therefore believe that more research is required into the development of pseudopoten- 
tials for explicitly correlated methods such as QMC. 
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